IFUM 545/FT November 1996 



The 2n-point renormalized coupling constants 
in the 3d Ising model: 
estimates by high temperature series to order (3 17 



P. Buter J and M. Comi0 

Istituto Nazionale di Fisica Nucleare 
Dipartimento di Fisica, Universitd di Milano 
Via Celoria 16, 20133 Milano, Italy 



Abstract 

Abstract: We compute the 2n-point renormalized coupling constants in the 
symmetric phase of the 3d Ising model on the sc lattice in terms of the high 
temperature expansions 0(/? 17 ) of the Fourier transformed 2n-point connected 
correlation functions at zero momentum. 

Our high temperature estimates of these quantities, which enter into the 
small field expansion of the effective potential for a 3d scalar field at the 
IR fixed point or, equivalently, in the critical equation of state of the 3d 
Ising model universality class, are compared with recent results obtained by 
renormalization group methods, strong coupling, stochastic simulations as 
well as previous high temperature expansions. 
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I. INTRODUCTION 



In recent times a considerable effort has been devoted to the evaluation of the 2n-point dimen- 
sionless renormalized coupling constants (RCC's) at zero momentum for the Ising model in three 
dimensions. These quantities are of interest for constructing the field theoretic effective potential 
]l]j2| of a 3d scalar field at the infrared fixed point or, in statistical mechanics language, for the 
formulation of the critical equation of state of the 3d Ising model universality class f3j-|| . The com- 
putational methods, which so far have been used, include various approximate forms |6-11] of the 



renormalization group (RG), the field theoretic strong coupling expansion M, the high temperature 
(HT) expansion |||4|,|l^-[l4| and (single- or multi-cluster) MonteCarlo techniques pj[5-17|. 



In this note we want to discuss how helpful in getting a first estimate of the RCC's in the 



symmetric phase, can be extensive HT expansion data published long ago [18] and so far only 
partially analyzed. Indeed expansions as double series in the HT variables v = tanh{(3) and 
r = exp{(3H), where (3 is the inverse temperature, are available for the Ising model free energy in 
a magnetic field H on various 2-, 3- and 4-dimensional lattices. In particular, in the 3d case the 
series extend up to order v 17 for the sc lattice, up to v ls for the bcc lattice and up to v 10 for the 
fee lattice. By computing the 2n-th derivative of the free energy with respect to the magnetic field 
at zero field we readily obtain the HT expansion of the Fourier transformed 2n-point connected 
correlation function at zero momentum (also called the 2n-th susceptibility) 

X2n{v) = ^2 < S(0)s(x2)s(x 3 )..s(x2n) >c (1) 

X2,X 3 ,..X2n 

These expansions together with that of the second moment correlation length 

f 2 (v) = J? 2 ^\ are the essential ingredients for the calculation of the RCC's. 

s v ' 6x2 (v) ° 

The expansion of the second moment of the 2-point correlation function fi^iy) on the sc lattice 



has been recently extended in Ref. l£j. 



In terms of these quantities, the first few RCC's, in the symmetric phase, are defined H as the 
values g^, (n > 2), that the following expressions 

9i{v) ~ 



4! e(v)xi(v) 



96{V) = -^-(- 76? \ 3/ N + 10 1%( \ 4/ J 



V 3 f Xs(v) | 56 X6(^)X4(^) 2gQ xl(v) 



»! V e(v)xl(v) e(v)xUv) e(v)xt(v) 



/a V 4 ( Xiofo) , 1on X8(^)X4(^) 10ft xl{v) 



iop e 2 {v)xi{v) e 2 {v)x» e 2 {v)x» 

-4620 li {V \ X l {V \ + 15400- Xt{v) ' 



e 2 (v)xl{v) i l2 (v)xl(vY 

take as v | v c . The volume V per lattice site has the value 1 for the sc lattice, 4/3\/3 for the 
bcc lattice and l/\/2 for the fee lattice. 
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We recall that scaling implies that, as the critical temperature is approached from above, we 
have X2n — B2 n (v c — ■y)~'>' _ ( 2n_2 ) A , where A is the gap exponent. If we also assume the validity of 
hyperscaling, we have 2 A = 3^+7 (where v and 7 are the critical exponents of £ and x respectively), 
so that the RCC's are finite (and universal) quantities. The quantities <?2n are expected |2(| to 
be of the form g2n{v) — g^n + ^2n( u c — v) d + ... as v | v c , where the dominant universal scaling 
correction exponent 9 has the value 6 = 0.50(2) [25|| for the 3d Ising model. 

By changing in the functions g2n{v) the variable v into y = ^ 2 {v), we obtain the strong coupling 
expansions, through the order y 17 , of the functions J2n{y) @ whose values at y = 00 give the RCC's. 

Let us add a few comments concerning the HT and the strong coupling series coefficients of 
the X2n on the sc lattice that we have tabulated, up to order v 17 , in the appendix together with 
the coefficients of the second moment of the correlation function ^{v) in order to provide the 
interested reader with all data we have used and thus make our calculations easily reproducible. 
The expansion for Xi was first computed |2l|] through v 17 using the data of Ref. fi~8fl , but only 
recently we could check it against a completely independent linked-cluster computation through 
the same order We should only draw attention to a minor misprint in the last 2 digits of the 
coefficient of v 12 as reported in Ref. [21]. Concerning the strong coupling expansions, we notice 
that in Ref. § 76 (y) has been tabulated, for any space dimension, through order y , while 7s (y) 
and 710 (y) through order y 7 only. A further significant extension of all these series can still be 
performed only then a complete check against an independent computation will be possible 
for the coefficients tabulated here. 

While this work was being completed, we became aware of a related work [14], also devoted to 
the analysis of the data of Ref. |l~8[| , and where also the low temperature side of the critical region 
is studied. We decided therefore to present only the part of our computation, mainly concerning 
the higher RCC's, which was not already covered by the very thorough discussion of Ref. [14]. 
In fact the availability of a longer HT expansion of £ 2 enables us to study also individual RCC's 
rather than only ratios among them, and moreover gives access to the strong coupling expansions. 



II. NUMERICAL RESULTS 

We shall now present our estimates of the first few RCC's as obtained from either the HT or 
the strong coupling expansions and discuss various "biased" or "unbiased" numerical procedures. 
In a first and straightforward approach we estimate y^n by evaluating at v = v c [22 24] 



near diagonal Pade approximants (PA's) of the quantity f2n(v) = g2n" 3 ( v ) which has a Tay- 
lor expansion in v. This procedure is not convenient for extrapolating gio(v), which changes its 
sign at some < vo < v c . In this case we should consider instead the expression {^) 6 gw(v), 
which also has a Taylor expansion in v. Thus, (biasing only v c ) we obtain the estimates: 
gt = 1.03(3), gt = 1.93(8), g+ = 1.53(36), g+ = -2.0(9). 

Here, as in the rest of this report, our estimates are given by a suitably weighted average over 
the results from the approximants using at least 14 series coefficients and the uncertainties are 
measured, conservatively, only on the basis of the spread of the results obtained from the highest 
approximants, always allowing also for the (much smaller) effects of the errors in v c and 9. 

It should be noticed that the central estimate of g\ obtained above is slightly, but significantly 
larger than the well established RG estimate g\ = 0.988(4) |25[| . 

This discrepancy leads us to investigate whether and to what extent these values are also 
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affected by a "systematic" error due to the non-analytic corrections to scaling which can spoil the 
convergence properties of the PA's. It has been suggested in Ref. @| that these corrections can 
be allowed for, or at least their effects can be significantly reduced, by performing the quadratic 
mapping v = v c (l — /|L~yLa ) with p = 2\/2 — 1. Essentially the same results are also obtained by 
using appropriately designed first order differential approximants [27] in which we can bias both 



v c and the scaling correction exponent 9. We arrive thus at our final set of estimates 

gf = 0.987(4), gf = 1.57(10), gf = 0.90(10), g+ = -0.71(35). (2) 

While the value of gf is only slightly lowered (and thereby closely reconciled with the most 
accurate RG estimates), the central estimates of the higher gf n are significantly altered and the 
uncertainties are reduced. Therefore it appears that our initial very simple numerical approach 
was rather inadequate and moreover we infer that the amplitudes Af n of the scaling correction 
terms increase with n. Finally, if we notice that the uncertainties of our estimates grow rapidly 
with the order of the RCC's, it is clear why, with the presently available series, we have to restrict 
our calculations to the gf n with n < 5. 

It is interesting also to study directly other quantities such as, for instance, appropriate ratios of 
the functions g2n{v) which do not depend on £ 2 and might be less sensitive to the scaling corrections, 
as a means to understand better the actual uncertainties of our numerical procedures. We have 
therefore considered the expression T± = g 4 6 ffa \ viv c and we have obtained the estimate Tf = 
1.75(5) neglecting the confluent singularity and, otherwise, Tf = 1.59(5). Analogously, we have 
also examined Tf = g 4 8 ffa \v\v c an d have estimated Tf = 1.29(43) by the first method and Tf = 

0.92(13) by the second, while for Tf = g?§£|i,T»c we obtain Tf = -0.7(7) and Tf = -0.35(20), 
respectively. All estimates of the T + are then completely consistent with the corresponding separate 
estimates of the gf . Notice that the Tf are simply related to the coefficients Fi of the small field 



expansion of the "reduced effective potential" computed in Ref. [10| as follows: Tf = 96-Fj 
Tf = 1728F 7 and T 3 + = ^^F 9 . 
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Let us also recall that long ago the sequence of universal amplitude combinations If r+ ^ = 



X2 ( v \xfr+4(v) \ v -\ Vc , r > 1, which are strictly related to the Tf , was introduced in Ref. [28] and, 



by using twelve term series [12], the first few if were estimated to be if = 7.73, if = 157.5, 
and Ig = 6180. (with no indication of error). Our estimates by using the direct PA procedure are 
1+ = 7.81(3), I? = 161.7(3), If = 6395.(21.), while if we allow for the scaling corrections, we 
find 

/ 5 + = 7.92(7), If = 165.(4.), I 9 + = 6809.(120.) (3) 

As it appears from the smaller difference between the results of the two kinds of numerical 
procedures, the and especially the if turn out to be less sensitive to the scaling corrections 
than the gf n and therefore we assume that they can be be determined with higher relative accuracy. 
It is therefore interesting to notice that from the above estimate of if we get the value gf = 1.62(6). 
Unfortunately however, at the present level of accuracy, the other simple relations among the higher 
Tf and the if, like Tf = ^{+lf - 567 f + 280) etc., which follow from the definitions of the gf n , 
cannot be used for improving the estimates of the higher Tf , and therefore of the corresponding 
gf n , by expressing them in terms of the If. For instance, Tf turns out to be a small difference 
between large numbers and the uncertainty of if is strongly amplified. For similar reasons it is 
also not useful to start directly with the critical amplitudes of the X2n- 



4 



An unbiased study of the RCC's can be performed starting with the strong coupling expansion. 
In Ref. H an elaborate extrapolation procedure was proposed which involves the dependence of 
the series coefficients on the space dimensionality. We have not yet computed this dependence up 
to order v 17 and therefore we cannot reproduce this procedure. We can, however, try the simplest 

2 

approach to evaluate 72n(oo), which consists in forming [N + 1/N] PA's to the quantity yi2n~ 3 (v) 
and in dividing them by y. This procedure is not very efficient and the only reasonably stable 
results obtained are: gj = 1.1(1), gt = 2.1(2), g£ = 1.9(2). 

We can also evaluate the ratios by diagonal PA's: again we find reasonable results only for 
Ti = 1.81(4). All these values are consistent with our previous first evaluation of these quantities. 

Alternatively, we can generalize a technique introduced in Ref. , which consists in inverting 

2 

the functions z<in = 72n 3 ™~ 3 (v) ( a ft er checking that the dependence of Z2 n on y is monotonic) and 
in determining g^ n from the value of z% n where y = y(z2n) diverges. This is conveniently done 
by forming PA's of the logarithmic derivative of y. The results are then: gf = 1.01(2), g§ = 
1.63(6), g^ = 1.05(15). As indicated above, these procedures cannot be used for computing gf . 

In conclusion, we believe that the general consistency among the results obtained by apply- 
ing suitable approximation procedures to various quantities with somewhat different properties 
corroborates our estimates in (|2|). 



III. A COMPARISON WITH OTHER ESTIMATES 

Let us now proceed to a comparison with the results already available in the literature. 

Our values in (|2|) for g\ and for g$ are not far from the estimates g\ = 1.018(1) and g§ = 
1.793(16) obtained in Ref. [jl4| from the analysis of the same HT series. A similar remark applies 
to the estimates g\ = 0.988(60) and g^ = 1.92(24) obtained in Ref. (jDJ from a sixteen term HT 
series. Our result for Jg" in (||) is also not far from the recent estimate Ie~ = 7.84(2) of Ref. [|l4}| . 

As to the strong coupling approach, we us recall that in Ref. Q the estimate g§ = 1.2(1) was 
obtained from an eleven term strong coupling series. 

It is also interesting to perform a comparison with the results obtained in the most extensive 



recent RG study [ 10 1 , by the fixed dimension (FD) expansion [25| up to five loop order, resummed 
by the Borel-Leroy technique combined with an appropriate conformal mapping. The estimate of 
gf agrees perfectly with ours, the central values <?g~ = 1.603(6) and g% = 0.82(8) agree with ours 
within ~ 2% and ~ 9%, respectively. On the other hand the larger disagreement about the value 
of gf should not be taken too seriously because, as noted above, the uncertainty which affects 
the calculation grows with the order of the RCC. Let us also return to a previous remark, in 
noticing that from the estimates of the Fi in Ref. [25| one arrives at the values I§ = 7.945(7), Ij = 



167.45(65) and Ig = 6718.(81.) , in very close agreement with our estimates in (||). Unfortunately, 
1^ and Ig are actually rather insensitive to the values of Fj and Fg. 

We also ought to recall that an independent calculation in the FD scheme gave the estimates 
g§ ~ 1.50 in the two loop approximation 0, g£ ~ 1.622 at three loop order with Pade-Borel 
resummation 11, and g§ ~ 1.596 at four loops p. On the other hand, from a 3 loop computation, 
values for g% have been obtained || which range from 0.68 to 2.71, depending on the resummation 
procedure. 

The approximate truncation of the RG flow equations studied in Ref. Q yields g\ = 1.2 and 
g§ = 2.25, which are both clearly larger than our values, although the ratio g§ jg\ =1.56 agrees 
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well with our estimate. An analogous, but lower order truncation of the RG flow equations [If] 
had given g§ = 2.40. 

The e = 4 — d expansion approach has not yet been pushed beyond order e 3 . ft has been used, 
in Ref. [1C], to produce the (rather large) estimates g± = f.f67, g§ = 2.30(5), g% = f .24(8) 
and gf = — f.97(f2). Notice however that, since also the estimate of g\ is unusually large, the 
corresponding values of the (or of the Fi) agree very closely with the results from the FD 
approach. The e expansion of was examined also in Ref. where by Pade-Borel resummation 
the estimate = 1.653 was obtained. 

Finally, we recall that the MonteCarlo simulations of Ref. || indicate g$ = 2.05(f5), which 
is not very far from our estimate, while the simulations described in Ref. [16 indicate the values 
g§ = 2.7(2) and g% = 4.3(6), significantly larger than both the RG results and ours. 

A summary of the present situation is presented in Table 1 which collects our estimates of the 
RCC's along with the corresponding ones obtained by other methods. 



IV. CONCLUSIONS 



We may conclude that, although the various computational approaches do not yet agree per- 
fectly, they do appear to converge to common estimates at least for the lowest RCC's. Therefore, 
in view of the difficulty of these calculations, we believe that the present residual discrepancies 
should not be overemphasized. The e expansion is certainly still too short, and perhaps, even for 
the FD expansions, a further extension would be welcome. The HT series presented here are not 
yet long enough, the more so the higher the order of the RCC considered. Indeed, we might argue 
that, at the order v s , the dominant contributions to the HT expansion of X2n(v) come from corre- 
lation functions of spins whose average relative distance is ~ s/2n, so that present HT expansions, 
in some sense, still describe a rather "small" system. Analogous problems of size also occur in 
stochastic simulations [6,16,17]. Therefore further effort would still be welcome to improve the 
reliability, the precision and, as a result, the consistency of the various approaches. 
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APPENDIX A: SERIES EXPANSIONS 

In the case of the sc lattice the HT expansion of the susceptibilities \2n ar e 

X2 (v) = 1 + 6v + 30v 2 + 150i> 3 + 726v i + 3510i> 5 + 16710i> 6 + 79494^ 7 + 375174k 8 + 1769686i> 9 
+8306862v 10 + 38975286D 11 + 182265822»; 12 + 852063558u 13 + 3973784886i> 14 
+18527532310»; 15 + 86228667894v 16 + 401225368086u 17 ... 
X4 (v) = -2 - 48v - 636v 2 - 6480u 3 - 56316v 4 - 441360»; 5 - 3208812v 6 - 22059120v 7 
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-145118844v 8 - 921726704v 9 - 5687262012v 10 - 34255147920V 11 - 202130397708v 12 
-1171902072144V 13 - 6691059944460v 14 - 37693869995312v 15 - 209838929195580v 16 



X6 (v) = 16 + 816v + 19920v 2 + 336720v 3 + 4518816v 4 + 51745680v 5 + 527187600v 6 



+4909918704v 7 + 42581232864v 8 + 348466330096v 9 + 2717492365392v 10 



+20347129869456V 11 + 147133138147872v 12 + 1032333377642448v 13 
+7054626581880336v 14 + 47100223055946160v 15 + 308027458769860704v 16 



Xs(v) = -272 - 23808v - 917376v 2 - 23013120v 3 - 437798496v 4 - 6852038400v 5 

-92654596992V 6 - 1117875129600v 7 - 12306018523104v 8 - 125633562017024v 9 
-1204105704419712v 1() - 10936791715557120v 11 - 94844317893543648v 12 
-789993027282411264V 13 - 6351007395937161600v 14 - 49478915100503151872v 15 
-374818460005448106720v 16 - 2768750733973561834752v 17 ... 
Xio(v) = +7936 + 1061376v + 59036160v 2 + 2049776640v 3 + 52252083456v 4 

+1067338759680v 5 + 18429925693440v 6 + 278749670360064v 7 + 3786553881275904v 8 
+47053476826003456v 9 + 542381843641961472v 10 + 5862580439606155776V 11 
+59934902216969609472v 12 + 583578982058859276288v 13 + 5442873762995091611136v 14 
+48857090955221240911360v 15 + 423771319439035687985664v 16 
+3563795335882672497655296v 17 ... 

The HT expansion of the second moment of the correlation function /j, 2 is 

M2 (v) = 6v + 72v 2 + 582v 3 + 4032v 4 + 25542v 5 + 153000v 6 + 880422v 7 + 4920576v 8 
+26879670v 9 + 144230088v 10 + 762587910V 11 + 3983525952v 12 + 20595680694v 13 
+105558845736v 14 + 536926539990v 15 + 2713148048256v 16 + 13630071574614v 17 .. 



The strong coupling expansions of the 72« to order y 17 are 



-3/2 



74 (y) = y — [l + 12y + 6/ + 48/ - 630/ + 7272/ - 83292/ + 957312/ - 11035662/ 



+127433528/ - 1472947908/° + 17036529504/ 1 - 197169806676/ 2 + 2283416559216/ 3 



-26463582511368/ 4 + 306946999598144/ 5 - 3563327123879550/ 6 



+41404188226284120/ 7 . 



76 ( y ) = -L- [l + 18y + 90/ + 48/ + 576/ - 8352/ + 114528/ - 1528416/ 



+19952712/ - 255983472/ + 3240722592/° - 40613845392/ 1 + 505052958336/ 



-6242882909472/ 3 + 76802505994224/ 4 - 941288338072752/' + 11501158664782008/ 



-140176233789711696/ 7 . 



— 9/2 

78 (j,) = ^ [l + 24y + 192/ + 576/ + 54/ + 6720/ - 113016/ + 1753632/ 

56 L 



-25771326/ + 364798032/ - 5028161232/° + 67958735808/ 1 - 904828659212/ 2 



+11905472505792/ 3 - 155154712361520/ 4 + 2006059450196288/ 



-25765180820314374/ 6 + 329050927608994224/ 7 . 



710(2/) = [l + 30y + 330/ + 1620/ + 3330/ - 1080/ + 67200/ - 1314720/ 



+22683000/ - 363847600/ + 5564033040/° - 82249187520/ 1 + 1185208196160/ 2 



-16740515134800/ 3 + 232658153938560/ 4 - 3190497478487440/ 



+43262377733737920/ 6 - 581022341984542560/ 7 . 
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TABLES 



TABLE I. A summary of the estimates of g^n by various methods. 



Method and Ref. 



JL 



_£s_ 



HT 



unpl. 

y 

is 



Strong 
HT 
HT 

RG FD-expans 
RG FD-expans. 
RG e-expans. 
RG approx. 
RG approx. 
Strong corml. 
MCM 

MC |ra 



0.987(4) 
1.01(1) 
1.019(6) 
0.988(60) 
0.987(2) 

1.167 
1.2 

0.986(10) 
0.97(2) 
1.02 



1.57(10) 
1.63(5) 
1.791(38) 
1.92(24) 
1.603(6) 
1.596 
2.30(5) 
2.25 
2.40 
1.2(1) 
2.05(15) 
2.7(2) 



0.90(10) 
1.05(9) 



0.83(8) 
0.68 - 2.71 
1.24(8) 



4.3(6) 



-0.71(35) 

-1.96(1.26) 
-1.97(12) 
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